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PREFACE 

This investigation uses the sensors of the Skylab Earth Resources 
Experiment Package to measure stratospheric aerosols. The data was acquired 
when the Skylab vehicle was in the Solar Inertial mode and the sensors were 
viewing the limb of the earth with the sun behind the sensors. 

The investigation was designed to produce altitude profiles of 
aerosol attenuation coefficients for several wavelengths. The analysis 
tecnniques were also developed to invert the data for particle size 
distributions and number densities; however, the actual data analysis did 
not reach that level due to unexpected complications discovered in the 
data which had to be treated. 

There were several contributors to this investigation beyond the 
Principal Investigator and the Co- Investigator. We acknowledge and 
appreciate the efforts of Nancy Pol ky who performed much of the actual 
implementation of the analysis techniques. We also appreciate the 
assistance of Dr. Bob Curran of NASA Goddard Space Flight Center who shared 
his software to read and calibrate the S191 data, and Dr. David E. Pitts 
who shared his software to calculate the effects of refraction. Ken Wahl in 
also assisted during the last few months of the investigation. We extend 
special appreciation to Bill Johnson of Lockheed Electronics Company who 
on several occasions provided us with information not contained in the 
supplied documentation. 
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1.0 INTRODUCTION 

The aerosol content of the atmosphere is related to the volcanic eruptions, 
industrial pollution, aircraft pollution, explosions on the surface of the 
earth, and perhaps other sources. Atmospheric scientists have long been 
interested in the distribution of aerosols in the atmosphere. This is re- 
flected in the vast literature that has grown up on the subject and in the 
considerable efforts that have been made to determine aerosol distribution 
including hundreds of aircraft, rocket, and balloon flights and many ground- 
based studies using searchlights and lasers. Some of this was motivated by 
the purely scientific objective of understanding our atmosphere. However, 
much of it was also motivated by a desire to make practical applications of 
this knowledge in astronomy, meteorology, and other fields. One of the more 
recent areas of interest is in remote sensing. Models based on the rather 
sparse data available give an aerosol optical depth of about 0.05 above 20 kni. 

An aerosol level four times this has been shown to cause error of several per- 
cent in the classification of ERTS data. It is quite likely that aerosol 
levels could vary by factors much larger than four. Thus there is considerable 
interest, both from a purely scientific standpoint and from a practical applica- 
tions standpoint in aerosol distributions in the atmosphere. 

The main objective of this experiment was to obtain quantitative measure- 
ments of the brightness of the earth's limb and to use these to study the 
distribution of aerosols in the atmosphere. 

Techniques were developed to use limb measurements (by 5190, SI 91 , and 
S192) to obtain aerosol attenuation coefficients, number densities and particle 
size distributions as a function of altitude. They are described in Section 
2.0. These methods were applied to SI 90, SI 91, and SI 92 data. Because of 
time and data quality limitations, it was decided to concentrate on obtaining 
aerosol attenuation coefficients rather than attempt the more complicated 
analysis required to obtain number densities and particle size distributions. 

The results are described in Section 4.0. 
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2.0 APPROACH 

2.1 ANALYSIS TECHNIQUES 

In analyzing the data the first task was to remove the noise insofar as 
possible. The techniques used to do this are described in Section 2.1.5. 

After this was done the data vvas analyzed using two general techniques: 

(1) comparing the measured limb brightness curves with simulated limb bright- 
ness curves calculated for various atmospheric models, and (2) inverting the 
data to obtain a vertical profile of the attenuation coefficients. The tech- 
niques for doing this are developed in this section. We also include in this 
section a description of a technique for inverting the data to obtain particle 
number densities and size distributions (Section 2.1.3). However, this 
technique was not used in the data analysis for the reasons given above. 

2.1.1 Simulated Limb Brightness Curves 

In this paragraph we derive a mathematical expression for the limb bright- 
ness curve that will allow us to simulate the brightness that would be observed 
by the Skylab sensors. This curve will, of course, depend on the aerosol 
distribution assumed for the atmosphere. All required data had a scattering 
angle of near 180°. The error caused by assuming the measurements to be 180° 
were expected to be small, and the task of modeling the slant path for nbn- 
180° scattering was complex by comparison. All models were therefore computed 

for 180° backscatter. We carry out calculations for a pure Rayleigh scattering 

12 

atmosphere and for two aerosol distributions found by Elterman. ’ Each point 
on a limb brightness curve corresponds to a line of sight through the atmos- 
phere (or sensor pointing vector) which traverses various levels in the 
atmosphere down to some minimum altitude x. Such a line of sight is the 
line O'C shown in Figure 1, where O' is the position of the sensor. Here we 
assume O' C is a straight line. In fact it is Slightly curved due to 
refraction in the atmosphere. This effect is considered in Section 2.1 ,6. In 
Figure 1, the point A is the point of closest approach to the earth's surface 
and point B is some other point on the Tine of sight, a distance +£ from A. 
Distances on the other side of A are negative. The altitude at point B is 
X f z. The radius of the earth is Rg.. This is the mean radius of the earth 
as defined in NASA documentation. 
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In order to compute the radiance at O' measured by a sensor pointed in 
the direction O'A we calculate the contribution of each point B along the 
line of sight and add all these contributions. Only first order scattering 
will be included and it is assumed that the incident solar radiation is 
parallel to O'C, i,e. , we consider the case of pure backscattering. Point B 
is illuminated by an incident solar irradiance of 

00 . 

-/ 6(u) du 

I e S, 

0 

where I is the solar irradiance at the top of the atmosphere and w(u) 

^ -1 
is the atmospheric volume extinction coefficient (units km ). The integral 

is along the path B O' and the upper limit is actually the distance AO'. How- 
ever, since B(j) is completely negligible (i.e., equal to zero in our models) 
for attitudes corresponding to the point 0' and above, we have replaced the 
upper limit by We shall assume that absorption is zero so 3 is equal to 
the scattering coefficient. 

The proportion of this incident irradiance that is scattered back towards 
the point O' is given by the value of the scattering diagram or phase function 
P(o) for 0 = For a more complete discussion of extinction coefficients 

and phase functions, see Chandrasekhar.'^ We shall write P or P{s., 180®) to 
show it is a function of U as well. The radiance I P(i?,, 0=180®) e"( 
scattered towards O' from the point B is attenuated on the path BO' by 
factor e“^ Integrating over all points B along the path we obtain 

the following expression for the radiance along this path 


B{X) = 1 P(£, 180®) 3(£) exp (-2 3{u) du) dZ (1) 

Equation 1 can be simplified if it is assumed that P(il, 180) is independent of 
position in the atmosphere (i .c., independent of Z and x). In that case one 
can easily show that 


IP, 

B(x) [1 - exp (-2 /“ 3(A) dA 


( 2 ) 
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In what follows v/e shall assume that g(A) depends only on altitude in the 
atmosphere. In that case = 3C^) so that 

r B(j^) dl = 2/“ 3(2.) d2,. (3) 

~oo 0 

From Figure 1 

2,^ = (R^ + X + z)^ - (R^ + x)^ = + 2(R^ + x)z. (4) 


(z + Rp + x) 

d2, = yr. : dZ (5) 

\/ z'^+2(R^+x)z 


Now 3(2.) is zero for all values of 2, except those corresponding to small 
values of X + z, i.e., values for which 


z « R^ 

X « Rg 


Therefore we can replace d2, in Equation (3) by 

dl p: 

V z 

The integral of Equation (3) can then be written 


r 3(2.) /“ 3(z + x) ^ /“ 3(z) 

O ^^0 y/z ^ ^ 1 


dz 


Thus we finally obtain the following expression for limb brightness: 




This is the expression that was used to calculate most of our simulated 
limb brightness curves. 


( 6 ) 


( 7 ) 


(?!) 
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2.1.2 Inversion of Limb Brightness Curves 

4 

The Weyl fractional integral of order 0 of 4> is defined by 
h(<f(f:), X, 0) = ^ (0 (f-xV^'^ dC 

It can be shown that the Weyl fractional integral and the Fourier transform 
are related by the expression 

Fe~^ ($(x); t} = t*^ {sgn t) exp (ir i ^/ 2 ) Fe“^ {h(^(C) tf 
where the inverse Fourier transform of f(x) is 

Fe~^ {f(x); t} f(x) dx = F(t) 


Then 


4»{x) = Fe t^ (sgn t) exp (■n:i(p/2) Fe"^ {h {${?), x,0);t} 


From Equation (8) 


itn 


I„P(x) 


I^P(x).2B(x) 


= Zjm r tIM dq 


Let 


F(x) 


r*" g(q) 


/ dq 


r(|-) X 


which is the Weyl fractional integral of order 1/2 of 3. Combining these 
last two expression yields 

1 


F(x) = -A Lr- in 


r(l) 2v^ 


I,P(x) 


IqP(xT-2B(^T 


3(x) = Fe {t^^^ (sgn t) exp (Tri/4) Fe"^ F(x);t] ; xl 


and 
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This expression serves to invert the limb brightness profile, B, which is 
contained in F(x) to find the altitude profile of attenuation coefficients. 
Note that the wavelength dependence of the coefficients is contained only in 
the brightness signature, therefore the resulting profile is for a single 
wavelength, and must be repeated to determine wavelength dependence. 

The function P(x) was found by evaluating Equation (3) with 
and without P(x) using values of 3 from a model of the stratosphere. A 
ratio of these results was then used to determine the function P(x) as it 
is removed from the integral (mean value theorem). It was found that the 
P(x) did not significantly vary with changes in the model. The error caused 
by the use of the function derived from the model on the measured data we 
therefore accepted as inconsequential. 

2.1.3 Inversion of Intensities for Number Density and Particle Size 
Distributions 

Scattering of radiation by aerosols which have a distribution of sizes 

can be represented by a Fredholm integral equation of the first kind assuming 

all of the particles have the same index of refraction. Using the Phillips- 

5 6 

Twoomey inversion method * the following development allows scattering by 
particles differing in index of refraction. Inverting the measured radiation 
from such a set of scatterers requires simultaneous integral equations for 
unknown size distributions with an independent measurement for each equation. 
Each size distribution then describes all of the particles having a common 
index of refraction. 

Let the single integral equation be represented by 

g(x) = K(x,y) f(y) dy (9) 

a 

where g(x) is the measured radiation, f(y) is the aerosol size distribution, 
and the kerneT is the scattering function. Constraints are then imposed on 
the error, e, inherent in the measurement of g(x), and the unknown function 
f(y), such as 

2 

Z 6. £ e, some chosen constant (10) 

i ' 
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(F‘ (y))^ dy = min {11 ) 

a 

where F'(y) is the second or third derivative of F(y). 

A trial solution, P(y), may also be used, 

(f(y) - P(y))^ dy = min (12) 

a 

Through some suitable quadrature the integrals become matrices, and an expres- 
sion is obtained which incorporates the constraints through Lagrangian 
multipliers. The problem then is to minimize the expression 

(G - AF)'’’ (G - AF) + Y, (BF)’’’ (BF) + (F-P)''' (F-P) 

T 

where A = wK, w being the quadrature weights and A is the transpose of A. 

The matrix B is determined by the numerical approximation to Equation (11 )• 
Matrix differentiation with respect to F yields 

(AF - G) + -y^ HF + (F-P) = 0 

where B = H, and solving for F yields the solution 

F = (a’*' A + H + y^)'^ (A*^ G + Y 2 P) 03) 

When two sets of particles which differ by index of refraction, each 
having an independent size distribution, are interspersed in a volume, the 
measured intensities scattered from the volume are due to the combined scat- 
tering of each set. If the intensities add linearly, which is the case for 
the scattering of light from an incoherent source, the problem of finding 
the unknown size distributions can be expressed in two simultaneous integral 
equations. 

’g(x) = /'’ ”K(x,y) '' f(y) dy ■. ^K(x.y) ^f{y) dy (14) 

a a 
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^g(x) = ;*’ ^''K(x,y) ’f(y) dy + /*’ ^^K(x,y) (y) dy 

a a 

N. . 

In these equations superscripts are used to represent different functions. 
Again the g(x) functions are the measured radiation* the kernels are the 
scattering functions and the f functions are the unknown particle size 
distributions differing in index of refraction. As before, the equations 
will be constrained through imposed limits on the inherent errors in the 
measured values and on the behavior of the unknowns in direct analogy to the 
solution of the single integral equatiori, A matrix expression is obtained, 




(^G - ^h V - (^G - ^’aV - ^F) + 


Y,(^bVi’'’ (’b’f) + Y 2 rB'=F)'(‘=B^F) + V 3 CF - 'P)' ('F - 'P)+Y 4 (‘F-‘P) ' (^F-‘P) 

1 2 

which is minimized with respect to the unknowns F and F yielding two 
simultaneous equations. 

”a 1 "a’f+’Vf- ’g) + 21 ;^T , 2242 p. 2 s).Y, 1 b^ 1 b 1 f^ 3 (’f- 1 p) = 0 


f2r>2r-\T/2n2 


Ir ln\r /It- li 


r2|- 2n ^T /2r- 2r 


^ V (^^A^F+^ Vf ~ ^G) + ^V(^^aV+^Vf^^G) + = 0 


These are then solved ^B) yielding 


/^ \ 


!■ 


llaT lU . 21.T 21 


n„T 12„ . 2UT 22, 


'AV ‘‘A +. 'A'_ 'A + T;i H + Y 3 I “A' '"A + "‘A' ■ "A 


\ ^2 / ^^^A^ ^^A + ^^A + ^^A + Y2 Y4J 


(15) 


^ ”a^ ’g V^V>G + T 3 I ^ 

\1V 1g + 2g + Y 4 I 
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The method is quickly generalized to n equations with n unknowns. The 
limitations are the same as those for solving the single equation with the 
additional requirement that measurements and hence the kernels be independent 
functions. The general expression to be minimized is 


N . 


?: j(^G - X r ^*'*A’'F)j + E L(’’bV)^ {^b''f)+u-(V-^P)^(V-’p)1 

' i=l i=l ^ i=r V ^ ' 


Differentiation with respect to each of the unknowns yields 

V ( I V) + Ytr V ‘^F + m.(^F-*^P) - = 0;K=1, 2,...l 


1=1 


or 


>: ^*^A + (y,^ I) 5^. 


j: ^G + p. *^P;f<=l,2,...N 

j ^ 


(16) 


where = V S. 


Comparison with Equation (13), the solution for the inversion of the single 
integral equation, shows that a similarity in the form of the solutions is 
maintained. Equation (15) becomes an example of the expansion of Equations (16) 
necessary to obtain the solutions* in this case two solutions from two simul- 
taneous equations. 


Respective number density is obtained from the solutions for size 
distributions by integrating over all sizes, 

N = /P(a)da 

Solutions can also be verified by calculating attenuation coefficients 
with the expression 


B(A) = /Cg^,g(A,a) f.(a)da where 

*^sca 3^ scattering cross section for a spherical particle of radius a. 
Again f(a) Is the particle size distribution. This expression can also 
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serve to invert the attenuation coefficients to obtain the size 
distribution. 

2.1.4 Location of the Field of View 

Locating the sensor field of view was a major task. The Skylab 
ephemeris and ancillary data did not give an accurate enough determination 
for our purposes. From NASA Document PH0-TR524 Rev. A, page B-1 , location 
of the intersection point of the sensor, pointing vector and the earth 
ellipsoid was accurate to within 5.5 Km, tolerance when the sensor was 
pointing vertically downward. Extending the tolerance from a distance of 
435 Km (orbital altitude) to 2400 Km (distance to the limb) gives a possible 
error of about 30 Km. Attempts were made to simply translate the field of 
view determined by Skybet, but all were unsatisfactory. Not only was the 
field of view unknown, but also the change in field of view between data 
points was unknown, as the following demonstrates: 

Let S^V be the sensor pointing vector in ECT‘ , a coordinate system 
very close to EOT, the Earth Centered True coordinate system, and let S^V 
be the sensor pointing vector in ECT. SC is the spacecraft vector in ECT. 

If A is a 3X3 matrix which accounts for the error In Sl^V, then 

A • sh' = sl»V 

If e is the angle between the spacecraft vector and the sensor pointing 
vector emanating from the spacecraft 

cos 0 = = (a ■ sfv)-tc 

[s?v| ■ |Scl ]S ■ sPv'[ |sc| 

defining h - spacecraft altitude 
R = radius of the earth 

X = altitude, or closest approach of the pointing vector. 
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X = (R+h) sin 0 - R 

= (R+h) sin {cos-' . R 

X 2 - X, (R+h) {s1n {cos’^ — j 

!a • SPV'I |scj 


- sin {cos ■ 


(t • SPV^) ‘ - SC 
- SPV'I Isti 


which shows that the difference in altitude is not a linear relationship 
with the difference in sensor pointing vectors. Even if SP^-j and SP ^2 are 
close the unknown transformation matrix % which contains the pointing errors 
cannot be removed. It is simple to show with the above relationship thi ' * 
I'', deviation from unity on the diagonal of A with very small {< 1%) off- 
diagonal contributions will produce the expected 30 Km error tolerance. 

It was for this reason that Skybet data was not suitable to use for 
fi el d-of-view location. 


2.1 .4.1 S192 Field of View 

We now derive an expression for the altitude of the sensor pointing 
vector for Sl'i2. (This expression was derived by Dr. F. E. Alzofon of 
Lockheed Electronics Co.) Each point on a SI 92 scan of the limb corresponds 
to a line of sight through the atmosphere (or sensor pointing vector) Which 
traverses various levels in the atmosphere down to a minimum altitude x. 

The situation Is shown in Figure 3. Here 0 is the center of the earth and 
Skylab is located at O'. The axis of S192 is pointed in a direction x' 
and the instrument is scanning in a cone as shown. The angular change 
from one data point to the next is Atj) 0.094 degree. O'P is a specific 
sensor pointing vector corresponding to a specific data point on a scan. 

In order to interpret the data we need to know the minimum altitude r 
measured from the center of the earth for this (and every other) sensor 
pointing vector. For the purpose of this discussion we will ignore the 
effects of refraction. They vnll be discussed in Section 2.1.6. 
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Figure 3. S192-"ExuCt Location of Field of View, A 


In order to calculate r we first define two coordinate systems centered 
at points 0 and O' with axes as shown. The axes in the two coordinate systems 
are parallel to each other. In the coordinate system centered at 0, the 
equation of the line O'P is given by 


X - x^ _ y ~ y-j _ ^ 

r™ ~ ~m ' ” n ~ 


(17) 


where x-| , y-j > and 2 ^ are the coordinates of the point O' and are given by 


^1 Vr^ - d^ 

y^-o 


(18) 


and f, m and n are the direction cosines of the line O'P. To determine 1, m 
and n, we consider Figure 4. Clearly 


REPRODUCIBILITY OP THE 
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Figure 4. S192--Exact Location of Field of View, B 

1 = cos a 
m = sin a sin (J> 
n = sin a cos ([> 

The equation of the sphere of radius r Is 

2 2 2 2 

Solving Equations (17) and (20) together we obtain 

A y^ t B y + C = 0 
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where A = 1 + a + b 
B = a •*- 2 db 

C = - r^ (22) 

a = 1 / tan a sin if) 
b = 1 / tan ^ 

Now Equation (21) has two solutions corresponding to the two points in which 
a line Intersects a sphere. However, we are concerned with the case where 
these two points are the same point, i.e., vsihere the Tine is tangent to the 
sphere- In this case 


- 4AC = 0 


Substituting Equations (22) In (23) we obtain 

r = - U-] a + d sin a cos if))‘ 

where x-j = - - d^ 


(23) 


The altitude of the sensor pointing vector is given by r - R^. 

A slightly different (and simpler) approximation to the above exact 
expression which does not depend on d (see Figure 3) can be obtained in the 
fol lowing manner. 

For S192 the altitude increment is dependent on the position of the scan 

relative to the earth. This can be represented as 
the angle of intersection, n-j a between the scan arc 
and the earth's surface. All altitude increments 
can be determined from the n-j a>id known charac- 
teristics of the scan arc. The arc between two 
adjacent points subtends an angle of a = 116.25°/ 
j 1239 pts. The arc between the first point and any 

Figure 5. SI 92— Altitude ^ 

Increment from Geometry pomt i subtends an angle of 
of Conical Scan 



0. - (i-1) (116.25°)/! 


239 
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To express in terms of n-j and 3^- it is seen in Figure 6 that n- = n^. -F , or 


”i ° "l {116.25®) 


Therefore the change in increment is 


AX. - S sin n. 

1 'i 



Figure 6. SI 9i!— Intersection of Conical Scan with the Earth's Surface 


where S is the center-to-center spacing between points i and i + 1. S is the 
product of the angle subtended by the arc between two adjacent points and the 
distance from the sensor. The angle (,157894 mill i radians) is a sensor 
characteristic; the distance from the sensor is a function of the altitude 
but almost a constant; a distance of 2400 Km was used. S then equals .378 
kilometers. For the cases where the arc reaches a maximum altitude and begins 
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decreasing in altitude, that Is,, at the maximuni, the tangent of the arc is 
parallel to the tangent to the earth; n| can be determined since the orienta- 
tion of the scan is known* For the cases where such a maximum cannot be 
identified in the data, n-j must be allowed to vary as the data is compared 
to the model thereby determining the unknown constants through a least- 
squares fit of the two curves. 

2.1.4 2 S191 Field of View 

The observation angle and the observation altitude for SI 91 ms derived 

to be 

0 = cos"V (sin Y cos 3) 

X = Rq sin 0 “ Rg 

as follows: 

The angle 0 is defined as the angle between the sun vector and its 
projection into the orbital plane. The angle y is defined as the sun 
elevation angle at the spacecraft subsatellite point. The angle y is also 
the angle between the sun vector and its projection into the orbit's tangent 
plane at the spacecraft. 

Some assumptions were necessary. The sun vector was not parallel to the 
EREP axis vector by a small angle less than four degrees. For this deriva- 
tion the two vectors were assumed parallel. Also the S191 vector is treated 
as equal to the EREP vector for this derivation of increment change between 
observation angles; the actual location was not the issue. Location was 
included in the derivation of scattering angle. 

Where 0 is the observation angle relative to the radius of orbit, 
from Figure 7, 

|DBI = Rq cos (90 - y} 

= Rg sin Y 
|ABj = jDBi/cos 3. 

eeproducibility of the 

OEIGMAL PAGE IS POOR 
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Therefore, 


lABl == Rq s1n y/cos B 
I AD I - ]ABl sin 3. 

Thev’efore, jADj = sin y sin 3/cos 3. 
|DCj - Rq sin (90 - y) 

= Rq cos y 

and jAC|^ - )AD|^ + ]DC!^ 



SI 91 — Location of Field of View 


then 

Also 


|AC|2 = r2 + r 2 cps2 T. 

° cos 3 

1AC|^ = iBAl^ + |BC!^ - 2|AB|[BCl cos 6 


where jBCj Rg‘ 

Solving for cos 0 and substituting the above relationships will yield 

cos 0 = sin Y cos 3» 

In this relationship 3 is nearly a constant and y is a function of time 
due to spacecraft orbital motion. Let L and M be the positions of the space- 
craft at times T-] and T 2 representing the acquisition times of two adjacent 
data points. The times are clearly related to the sensor scan rate by 

^2 " ^1 ^ ^scan rate^ 

= T-j + . 925 sec 

The angle, a, defined by the extended earth radii passing through L and M 
respectively, is determined by the product of the orbital rate and the 
inverse of the scan rate. 
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If point N is defined by the intersection of their respective orbit 
tangent planes vjith the orbital plane and S^M and S^L are parallel 
projections of sun vectors into the orbital plane, then 

Yg "" Tl " “ 

The initial observation angle 6 ^ is related to the initial estimate of 
lowest observed altitude, 

e, = sin''' [(Re + 

As previously shown 6 and y are related by 



sin Y cos 3 = cos G 

then 

Y-| = sin (cos 0 ^/cos 3 ) 

and 

7 i+i=^i'« 

and 

0^. = cos"^ (sin Y- cos a) 

and 

X. = sin 0 . - 

where 

is the altitude of point i. 
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2. 1.4. 3 S190A Field of View 

The field of view of S190A as described in NASA PH0-TR524 is 21.24° x 
21.24° which is about 900 Km at the earth's limb 2400 Km distant 
(2 J2400j tan ). The data frame Is 57 x 57 mm. Using a micro- 

densitometer with a sensing aperture of 20 y, each data point records .316 

Kjn { - £ QP , at the earth’s limb. If data points are recorded 

57 X 10'’ y 

every 20p, center to center spacing is also .316 Km. 

The S190A densitometer scans were aligned to be approximately normal to 
the earth. If the alignment were in error up to 20°, the altitude increment 
would be between .30 Km and .32 Km. 

2.1.5 Noise Reduction 

Noise elements identified and treated consisted of a Gaussian or random 
noise, a high frequency noise, a low frequency noise, and a coherent noise 
generally known as "herringbone" noise. 

The random noise was to be expected. Nevertheless the random nature of 
the noise was verified by constructing noise histograms and verifying their 
characteristic Gaussian shapes. The histograms were not exactly Gaussian 
due to the contributions of coherent noise and aerosol layer effects, however, 
the curves approximated the Gaussian shape well. 

The random noise was obtained for histogramming by fitting the measure- 
ments on the limb with a sixth order polynomial and finding the deviations 
from the polynomial. This removed the effects of amplitude. Using a 
polynomial approximation, however, does not account for the "bumps" caused 
by aerosol layers, therefore the Gaussian shape of the histograms should 
have a small distortion as they did. 

The random noise was reduced by averaging twenty- five adjacent scans. 

The number of scans to average was arrived at by considering the advantage 
of using more scans (the standard deviation decreases as the inverse of the 
square root of the number of scans) versus the rate of field of view shift 
across pixels. The shift rate was found by determining the average apparent 
point of earth-atmosphere Interface. Two such points were found 2692 scans 
apart. With a change in the pixel number of 300, the rate becomes about 
.11 pixels per scan, or 2.8 pixels per 25 scans. 
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A low pass Fourier filter was applied to the result of averaging. The 
filter was designed to suppress high frequency noise equivalent to a period 
on the limb of less than about 2 Km. The filter was applied to deviations 
of the data (average curve) from a sixth order polynomial least squares 
approximation to the data. Figure 9 shows the typical result; the 
"herringbone" noise reported by NASA becomes obvious. The frequency 
measured in this data appears to be about 10 KHz. The Fourier filter was 
therefore modified to suppress this noise, which reduces the effective 
resolution to about 2.2 Km on the limb. 

After the data had been Fourier filtered it was apparent that we were 
not successful in removing all coherent noise. A twenty-five point convolu- 
tion filter was therefore applied to the data in the time domain rather than 
Fourier transform the data again to be filtered In the frequency domain. 

This filter removed much of the noise; however, some lower frequency 
noise remained. For this reason the sixth order polynomial fit of the data 
was analyzed in addition to the filtered data; this givev a general level of 
aerosol content. 

The low frequency noise was suspected when it was noted that the least 
squares sixth order polynomial approximation to the data showed a shift in 
amplitude between several adjacent scans. A 17 cycles/ sec noise was re- 
ported by NASA personnel. With a scan rate of 94.8 scans/sec we searched 
for a noise with a period of approximately 5.5 scans/cycle. 

The frequency and amplitude v;ere determined by plotting a selected point 
on a scan as a function of scan. A twenty point curve was thus derived, for 
example the value of data point 500 could be plotted for scans 200 through 
220. Several such curves were obtained by repeating the procedure for 
different points on the scan. To emphasize differences and to minimize 
effects of amplitude, the lowest value on each twenty point curve was sub- 
tracted from the respective curves. The plots for each set of scans were then 
averaged to obtain a representative curve for that set. Plots were obtained 
both from the portion of the scan covering the limb and from the portion 
of the scan record!, j dark space. The limb plots are important simply be- 
cause that is the data being analyzed for aerosol effects; it is preferable 
not to extrapolate results from other data into this data. The dark space 
plots are valuable since Ideally they record no external energy; any 
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Figure 9. S192--Pass 47, Tape 934527, Scans 350-374, Band 3. 

Deviations from Polynomial Fit After Fourier Filtering 
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fluctuation from scan to scan should be due to the sensor alone thereby 
making the variations more apparent. For the limb plots the value plotted 
was obtained from the polynomial approximation of the data; for the dark 
space plots the value plotted was obtained from a straight line approximation 
of the data. This was necessary since the Gaussian noise would camouflage 
the periodic noise being sought. 

Figures 10 through 12 are representative of the low frequency plots for 
the limb. Figure 13 is representative of the plots derived from dark space. 

A comparison shows that they are consistent both in amplitude and in 
frequency. Both sets reveal a frequency of about 5.8 scans/cycle or 16.3 
cycles/sec. 

The low frequency noise was therefore determined. The result on each 
scan is an additive constant which is unknown, but related scan to scan. 
However, the above procedure to suppress random noise averages twenty-five 
adjacent scans which covers about four cycles. Therefore, only an unknown, 
additive constant remains. This constant is determined by a solution of 
an overdetermined set of equations during the scaling process. A more 
complete description of the scaling process is found in Section 3.1.2. 

2.1.6 Refraction 

Refraction was included in the scaling of the data by assigning the 
appropriate observation altitude to each measured point. The model remained 
that which would be observed without refraction; the location of data points 
were distributed to remove the effects of refraction. 

To develop the proper relationships the increment change in observation 
angle on the vertical was first determined; this increment is not that incre- 
ment inherent in the design of the instrument as described in Section 2.1.4. 
Since we are seeking to describe an equivalent sensor which observes a path 
normal to the earth one need not consider what is observed but simply the 
observation charaeteri sties. The geometric relationship of Figure 14 yields 

0^ = sin"'' [(Rg + X^)/R^] 
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Scan 

Figure 10. S192— Pass 61, Tape 932867, Barifi'*!. Low Frequency 

Scan to Scan Noise for Four Points on the Limb 



Figure 11. S192--Pass 61, Tape 932867, Band 1. Low Frequency 
Scan to Scan Noise for Fo,ur Points on the Limb 
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Tape 932857 


SI 92— Pass 61, Tape 932857, Band 3. Low Frequency 
Scan to Scan Noise for the Average of Deep Space 






Research fi Engineering Division ^ 

BOEING AEROSPACE COMPANY JT^ 



Figure 14. Observation Angle at the Sensor for a Straight Line Path 

Thus for SI 92, an altitude provided by the derivation of Section 2. 1.4. SI 92 
can be converted to an equivalent observation angle 6. For ST91 Section 
2.1.4.S191 derives the observation angle directly. Knowing the observation 
angle, one can obtain the altitude viev;ed by tracing the refracted light 
path through the atmosphere. The data point acquired at the given observa- 
tion angle is then assigned the obtained altitude, the complete set of 
which defines the limb brightness profile with the effects of refraction 
removed. This profile is then compared to the model for scaling and calibra- 
tion of the data as defined in Section 3.1.2. 

To trace the refracted light path the atmospheric refraction model as 
implemented by Dr. David E. Pitts of NASA and documented in NASA TMX-58033, 

7 

"A Model Atmosphere for Earth Resources Applications" was used. For a 
path leaving altitude X at an angle y from the earth's normal, this program 
uses the "law of sines" to find 4>*, the angle between the path and the 
earth's normal at an altitude Z above X. The program considers the atmos- 
phere as a number of layers and iterates through them until the light path 
emerges from the final layer. The geometry is illustrated in Figure 15, the 
relationship for refraction is 



The indices of refraction, n' and n", of the respective layers were computed 
within the same program. Because refraction was only a small effect, any 
inaccuracies caused by approximated input parameters would be insignificant. 
Therefore, one set of temperature and pressure as a function of altitude 

Q 

were used throughout, being provided by the U.S. Standard Atmosphere, 1962. 
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Because this investigation observes the lowest altitude at a tangent 
(> = it/2) a table of altitude versus emerging angle was derived by assigning 
an altitude and initiating Y to be tt/2 and iterating through the atmosphere 
to obtain the emerging angle. Because essentially no refraction occurs above 
20 Km iteration was terminated there. The emerging angle, 4>', was related 
to the observation angle, 0, at the spacecraft by the expression 

8 = sin"^ |(Rg + 20) sin (|)7 Rq| 


as can be seen by Figure 16 is the "law of sines". 

The program was verified by using another expression derived by Baum 
g 

and Code which derives the amount of refraction occuring as light passes 
through a planetary atmosphere having exponentially decreasing density as 
altitude increases: 


♦l ' - l) [(R' - Rp/H 
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where n' is the refractive index at X', a reference altitude for which n‘ 
known (R‘ = Rg + X'). X- - R' - Rg is the altitude of closest approach 
to the planet; is again the radius of the earth. A scale height, H, of 
8 Km was used. 

The Pitts program calculates (}j' and a' , the Baum and Code expression 
calculates shown in Figure 17. It can be seen that 

+ ^/2 - a. 

Comparison of the results of the two agree to better than "it which is as 
accurate as ti.a value used for the scale height parameter. 
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2.2 DATA ACQUISITION 

This investigation used data from S190A, S191, and S192. All data was 
acquired with the Skylab vehicle in the solar inertial mode as opposed to the 
Z-Local “Vertical mode. This positions the spacecraft such that the sensors' 
pointing vector views the earth's limb twice per orbit. Each of these two 
limb passes provides a complete limb data set. 

2.2.1 The S190A Data 

The S190A sensor consists of six 70 mm Innses mounted in a single 
camera body. It produces film Imagery of a scene in six spectral regions. 

The regions and film types are listed in Table 1. 

Table 1. 190A Spectral Bands 


Spectral Region 

Film Type 

,5 - .6 pm 

B/W SO-022 

.6 - .7 

B/W SO-022 

.7 - .8 

IR B/W EK2424 

.8 - .9 

IR B/W EK2424 

.5 - .88 

IR Color EK2443 

.4 - .7 

High Resolution Color 


The S190A camera bank was used because it provided maximum spatial 
resolution and minimum geometric distortion of the image. 

2.2.2 The S191 Data 

The SI 91 spectrometer was used to acquire high spectral resolution 
data within its field of view of one mllliradian (about 2A in altitude 
at the earth's limb). Radiance data from the first filter wheel was used. 
This data includes wavelengths from 0.38 pm to 2.5 pm and has a wavelength- 
dependent spectral resolution of approximately 0.08 X. We did not use the 
SI 91 data in the infrared bands and the first, or lowest, wavelength bands 
since it was reported that this data had errors dfj.i to off-axis rejection 
problem. An algorithm was developed to correct for this effect but we were 
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unable to acquire data corrected by this algorithm. 

Approximately one spectral scan was acquired per second. In solar 
inertial mode the orbital motion of the satellite caused the line of sight 
of SI 91 to traverse the altitude at the earth's limb. The data was there- 
fore acquired by fixing the sensor pointing angle relative to the spacecraft 
and allowing the spacecraft motion to produce the altitude scan. At one 
spectral scan per second the line of sight traversed about 2.4 Km in 
altitude when the motion of the sensor line of sight was normal to the 
earth and in the orbital plane. For a line of sight off of the orbital 
plane at an angle of about 44°, it traversed about 1.7 Km per spectral scan. 

The SI 91 spectrometer recorded three complete scans through the limb 
for each complete limb data set. This was to obtain data in the same ap- 
proximate region of the earth for three angles of limb observation relative 
to the direction of the Incident solar radiation, all near 180° backscatter. 
The sensor was pointed just below the horizon of the spacecraft to cause 
the field of view to scan upward through the atmosphere. At intervals of 
(35/cos 3) seconds the VTS'^was manually reset to below the horizon and the 
process was repeated. The angle, 3> is that angle defined by the inter- 
section of the sun to earth center line and the projection of this line on the 
orbital plane. The first scan was Initiated when the spacecraft was about 
33° from the terminator on the daylight side and repeated until the vehicle 
was about 13° from the terminator. These guidelines were derived from the 
fact that near 33° from the terminator the limb becomes visible with the VTS 
about 10° forward; at 13° from the terminator the limb is seen behind the 
spacecraft at angles greater than 10°, It is during this period, at a time 
where the spacecraft position in orbit is 3.25° before the terminator, 
that the S192 scanner initiated its data acquisition for approximately 
five seconds. At least three frames of S190A were exposed during this 
period al so. 

2.2.3 The SI 92 Data 

The third sensor of EREP used in this investigation was the multi- 
spectral conical scanner S192. Data was recorded simultaneous Ty in thirteen 
spectral bands with wavelength between .41 p and 12.5 y as shown in Table 2. 
The conical scan has an angle of revolution, a, of 5°32‘ . Experiment 
data is acquired in the front 116° 15’, the remaining portion of the scan 

^Viewfinder Tracking System • 

■ - SI ■ . ■■ 
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Table 2. SI 92 Spectral Bands 

Band 

1 

2 

3 

4 

_ 5 

6 

7 

8 

9 

10 

11 

12 
,13 

being assigned to housekeeping data. One conical scan was recorded approxi- 
mately every .01 second, therefore the effect produced by motion of the 
spacecraft during one scan was negligible. The scanner was allowed to 
record several thousand limb scans as it passed through the limb. This data 
provided a horizontal distribution of stratospheric characteristics as well 
as sufficient data to reduce system noise effects. 

2.3 DESCRIPTION OF ACQUIRED DATA 

Analysis of the S190A data was done for Frame 175 of Pass 47 taken on 
September 18, 1973 at 44 minutes 8 seconds after midnight, GMT. Three micro- 
densitometer scans were taken from the frame. The first was taken one- 
half inch from the left edge of the frame, the second was taken in the mid- 
dle of the frame, and the third was measured one-half inch from the right 
edge of the frame. Table 3 lists the variable information for these three 
scans. 

Only one 5191 scan was analyzed. This was the first data take of Pass 
47 taken on September 18, 1973. The scan covered from 44 minutes 19 seconds 
through 44 minutes 52 seconds after midnight, GMT. The sun elevation 
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angle at subsatellite point on earth was 2.63^ and the data was acquired 
at 43°N, 104^J. 

Eight sets of scans were analyzed from the 5192 data. Table 4 sum- 
marizes the variable information of these scans. The five sets of scans 
from Pass 47 were taken on September 18, 1973, 44 minutes after midnight, 
GMT. 

The three sets of scans from Pass 61 vnere taken on December 5, 1973, 
16 hours and 45 minutes after midnight. 


Table 3. S190A - Pass 47, Frame 175 



Latitude 

Longitude 

Sun-angle 

Left 

47° N 

103° W 

2.31° 

Middle 

45° W 

103° W 

2,21° 

Ri ght 

42° N 

103° W 

2.16° 


Table 4, 

SI 92 



Pass 

Tape No. 

Scans 

Time (sec)* 

Lat-Long** 

Sun-angle 

47 

934527 

350-374 

10.80 

43°N-104°W 

2.58° 

47 

934527 

600-624 

13.44 

43°N-104°W 

2.49° 

47 

934528 

100-124 

15.67 

43°N-104°W 

2.21° 

47 

934528 

300-324 

17.78 

43°N-104°W 

2.09° 

47 

934528 

614-638 

21.10 

42°N-103°W 

2.00° 

61 

932867 

100-124 

41.12 

45°S-33°E 

4.49° 

61 

932867 

600-624 

46.44 

45°S-34°E 

4.23° 

61 

932857 

120-144 

69.54 

45°S-36°E 

3.08° 


^number of seconds after stated times the scans took place 
**position on the earth where the earth's tangent Is parallel to 
the sensor line of sight at its closest approach (point a of 
figure 14) 
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3,0 DATA ANALYSIS 

3.1 THE S192 DATA 

The first step in analyzing the SI 92 data was to remove the noise 
insofar as possible using the techniques described in Section 2.1.5. This 
involved the following steps: 

(1) Put data in engineering units using the calibration constants 
supplied on the data tape lisader record. 

(2) Average 25 adjacent scans. 

(3) Fit the exponential of a 6th order polynomial to the result 
of step (2). 

(4) Find the difference between the results of (2) and (3). 

(5) Apply a Fourier filter to the results of (4). 

(6) Apply a convolution filter to the results of (5). 

(7) Add the results of (6) to the fit found in (3) for a smoothed 
data scan. 

In some of the SI 92 scans the maximum altitude was low enough that one 
could determine at what point in the scan it occurred. In other cases the 
scan extended so far into deep space that this could not be done. These 
two cases are treated somewhat differently. 

3.1.1 Position in Scan of Maximum Altitude Known 

In this case a relative altitude scale can be obtained using the 
methods of Section 2.1 .4. However, an absolute scale cannot be directly 
determined because the absolute value for the maximum altitude is not known. 

We attempted to determine the absolute scale by comparing the 
measured curve with a simulated curve calculated from Equation (8) of 
Section 2.1.1. By varying the altitude of the first point a least squares 
fit between the measured curve and the model curve was found as a function 
of starting altitude. The fit was found only in the last half of the 
curves, i.e., the higher altitudes, where the aerosol contribution is 
minimal and the measured data should be very close to the simulated data. 
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In such a comparison of curves, one expects to require some adjustment 
of the brightness as well as the altitude because there are uncertainties 
in both the measured data and the model. We assumed variations of the 
form = a T^- + b where the are the model values and the are the 
measured data values. For each choice of starting altitude, a and b 
are obtained through a least squares solution of an overdetermined set of 
equations. Scan points used for S^* and were chosen to weight the second 
half of the curves. Only every tenth point was included from the first half 
of the curves. Every point was included from the second half. It was hoped 
that a would be very close to one, and b close to zero. In that case, . 
one would have some confidence in the absolute altitudes obtained. 

The inversion of the data ms carried out in two ways. The first was 
to simply invert the measured data using the information of Section 2.1.2. 
Only the relative altitudes, which were known, are required for this. This 
leads to an altitude profile for 3 whose relative scale is known but whose 
absolute position in space is not known. This can be determined approxi- 
mately by relating it to the absolute altitude scale determined by the 
procedure described above or by comparing this profile to the one in the 
£1 terman model , 

The second data inversion technique involved inverting the data after 
it had been transformed using the coefficients a and b determined above. 

The process was the same except that the data had been scaled with a and b. 

3.1.2 Position in Scan of Maximum Altitude Not Known 

In the case where no peak is present, the analysis is complicated by 
the fact that there is no direct method for determining even the relative 
altitudes. 

In this case the data curve was also fit to the model by finding the 
least squares solution to the overdetermined set of linear equations, 
a ty+ b = S.j, where t^. is the measured brightness and is the model 
brightness at the same altitude. 

The altitudes were assigned using an iterative process which varies 
both with the assignment of an altitude to the first data point and with the 
assignment of the increment between the first two altitudes. As shown 
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earlier the first increment determines all subsequent intervals. An initial 
estimate of altitude was made based on the point of maximum slope of the 
curve. The point of maximum slope of the data is assumed to be at the same 
altitude as the model. This assumption could be in error by 3 Km but it 
serves well as a first estimate. By allowing the data to converge on a 
best fit the unknown constants are determined. An obvious shortcoming of 
this technique is the severe dependency on the model. In an attempt to 

1 o 

bracket the results the two Elterman ’ models of 1964 and 1968 were used as, 
respectively, a low and a high aerosol conten;: stratosphere. 

3.2 THE SI 91 DATA 

The first step in analyzing the S191 data was to average several 
adjacent intensity values within a spectral scan. Although this decreases 
spectral resolution it was necessary in order to reduce the random noise 
within each data value caused by low Irradiance. Table 5 gives the 
spectral bands created in this way 

Relative altitude increments were determined as described in 
Section 2.1.4; the absolute altitudes were not known. In the same manner 
as was described in Section 3.1.1 for S192, the data was compared to a 
simulated curve calculated from Equation (8) of Section 2.1.1 in order to 
determine absolute altitudes. Again it was hoped that the scaling constants 
determined by the least squares solution of the comparison would only 
slightly modify the calibration provided with the data by NASA. 

Table s, SI 91 Spectral Bands 


Band 

Spectral 

Range (ym) 

Average Wavelength (urn) 

1 

.46 

- 

.50 


2 

.51 

- 

.55 

.53 

3 

. 56 

- 

.60 

.58 

4 

.61 

- 

.65 

•63 

5 

.66 

*- 

.70 

.68 

6 

.71 

- 

.75 

.73 

7 

.76 

- 

.80 

.78 

8 

.81 

- 

.85 

,83 

9 

.86 

- 

.90 

.88 

10 

.915 

- 

.945 

.93 

11 

.96 

- 

.99 

.975 
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3.3 THE S190A DATA 

Calibration of film densities is done by constructing a density vs. log 
exposure curve from a calibration step wedge. These curves are given for 
each roll of original film. We have duplicate film, so certain steps must 
be taken to construct this curve. 

Machine calibration step wedges and the density measured at each step 
by PTD*=are supplied. There is a separate step wedge for each of the three 
duplicate film types. Using the densities that we measured on these step 
wedges, and those supplied, a machine cross-cal ibration curve was con- 
structed for each duplicate film type. With these curves, densities 
measured by PTD can be converted to equivalent densities which we would 
have measured with our densitometer. 

We were supplied with the densities measured by PTD on the step wedge 
for each duplicate roll of film. Using the machine cross-calibration curve, 
these densities were converted to the densities which we would have 
measured. Because each step on a duplicate step wedge corresponds to the 
same log exposure increment as the same step on the original step wedge, 
log exposure values for the densities we measure can be found by a user 
duplicate density vs. original log exposure curve. 

The original log exposure at each step is found from the log exposure 
at the first step of the wedge plus the log exposure increment for each 
step. The log exposure at the first step is supplied in the Sensi tometric 
Data Package. These values must be adjusted for conversion from the 
Wi^atten filter, used in the construction of the step wedges to the flight 
type filter. These adjustment factors are also supplied in the 
Sensitometric Data Package. The log exposure increments were supplied by 
Harold Lockwood of Technicolor Graphic Services, Inc. 

Log exposure is first converted to exposure, then exposure is converted 
to radiance values by the following equation 

D - ^ • 4 - 10~^ 

^ " f^Tt 

where 

R = radiance incident at the lens in Watt/cm 

E = exposure in ergs/cm^ 

*Photographic Technology Division - NASA JSC 
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F = f“Stop number 
T = optical system transmission 
t = duration of exposure in seconds 

Sensitivity is lost in the film duplication process. Figures 18 through 21 
show the original and duplicate characteristic curves for Rolls 43, 44, 47 
and 48 (Stations 1, 2, 5 and 6) of Pass 47. The loss of sensitivity is 
readily observable from these curves. As many as 10 density levels {steps 
of the step wedge) with a change in density of 1 on the original film 
produce constant density on the duplicate film. Only 7 steps of the step 
wedge or 'v^l.05 log exposure units lie within the linear portion of the curve. 

The S190A data was measured as duplicate film density with a 400 x 20 gm 
aperture densitometer oriented with length tangent to the earth. 

In addition to the convolution filter inherent in using a 400 x 20 urn 
aperture, the values were further filtered with a fourteen point convolution 
filter. This was necessary because the film grain noise still dominated 
the density fluctuations. 

As with the other sensors, the absolute altitudes and altitude increment 
were unknown. In the manner described in Section 3.1.1, the data was best 
fit to a model brightness curve by iterating on starting altitude and 
increment, to determine altitudes, scaled brightness values and attenuation 
coefficient profiles. 
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log exposure . . , 

^ original film 

duplicate film 

Figure 18 S190A--Pass 47, D-log E Curves for Original and 
Duplicate Film. Roll 43-24. 


Figure 19..S190A--Pass 47, D-log E Curves for Original and 
Duplicate Film. Roll 44-22 


log exposure 
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original film 
duplicate film 


log exposure 


Figure 21 . S190A-Pass 47, D-log E Curves for Original 
Duplicate Film. Roll 48-32. ^ 


Figure 20 S190A--Pass 47, D-log Curves for Original and 
Duplicate Film. Roll 47-32. 
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4.0 RESULTS 

SI 92 Results 

One of the more interesting results v/as observed in the SI 92 data 
(pages 49 through 82). A peak was apparent in the lower altitudes that 
was not present in the shorter wavelengths and grew with increasing wave- 
length beginning with band 7 (.825 pm). Figures on pages 49 through 72 

2 

are the radiance values in milliwatts/cm -pm-str from pass 61. Each has 
been averaged over twenty-five scans to remove most of the random noise. 
However, It is apparent that some random noise still remains. The first 
twelve charts are from scans 100-124; the second twelve charts aie from 
scans 600-624 showing the continuance for several kilometers. Pages 73 
through 82 show a similar peak from pass 47. These ten charts were scaled 
to an inaccurate model which therefore nullifies the radiance and altitude 
values shown. However, the charts serve to illustrate the similar wave- 
length dependent peak. Judging from the altitude values obtained from the 
scaling process previously described (and illustrated below) one can 
estimate the altitude of the peaks to be around 11 Km. There were cirrus 
clouds present in the S190A imagery. We therefore suggest that the peak 
could be a cirrus cloud at about 36,000 ft. The wavelength dependence is 
due to the limb view. The blue portions of the incident light are scattered 
out of the incident path. The cloud scatters red light which is further 
reddened through its path out of the atmosphere. Hence the sensor detects 
only the longer wavelength radiation. The shorter wavelength radiation 
effectively does not reach the lower level at which the cloud is seen. 

The charts of pages 83 through 87 are the attenuation coefficients for 
j ten SI 92 wavelengths obtained by inverting data from pass 61. The relative 

altitude increment was determined by knowledge of the relative position of 
! the highest point In the scan arc as was described earlier. The data was 

’ not sealed in intensity or altitude. The resul ts are therefore valuable 

! CO identify relative values and fluctuations in the attenuation coefficients, 

i Comparison vnth various model s or previous results would be necessary to as- 

sign altitudes, however assignment can be done by comparison with the fol- 
I lowing charts. To scale them accurately in this way would require an 

. additional least squares curve matching which was not performed. A 

I general interpretation is also possible simply by noting the attenuation 
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coefficient values and the relative variations. These variations may be 
due to sensor noise or to aerosol variations. We were unable to complete 
a study of these possibilities due to other complications. 

The charts of pages 88 through ‘159 contain the results of scaling 
and inverting the S192 data for passes 47 and 61 using the variable altitude 
and increment scheme as previously discussed. The data is presented 
by pass, for the two models. Each result is represented as three charts: 

(1 ) the limb brightness measurement as scaled with the model used to 
scale, (2) the attenuation coefficients obtained by inverting the brightness 
v/ith the model for a Rayleigh atmosphere (no aerosols) and with the model 
attenuation coefficients (the first and last 5 points should be disregarded 
due to deviations inherent in the Fourier inversion process), and (3) a ratio 
of the aerosol and Rayleigh coefficients to accentuate layers. 

Table 7 shows the additive and multiplicative scale factors applied 
to the data. Since one would expect the multiplicative factor, a, to be 
near unity if the NASA provided calibration were close, only those results 
which used factors close to unity are presented here. One would also expect 
that the additive bias, b, would be small compared to the data. The actual 
criteria used to select the best results were as follows: 

1) .95 < a < 1.1 

2) b < 0 and .5 < 

3) a good agreement in brightness values between the model and the 
measurements for the upper half of the curve. 

Among things apparent in the data is a layer centered at 40 Km about 
10 Km deep, a layer at 35 Km of approximately the same depth and one at from 
20 to 25 tdii from 5 to 10 Km deep. 

The 40 Km layer is apparent in the plots on pages 90, 99, 102, 111, 120, 
135 and 156. Of these scans, all but those on pages 102 and 120 meet the 
criteria established above. The plot on 120 meets the criteria except that 
jbj/deep space value = .35, slightly smaller than allowed. All of these 
plots are for bands 1 or 2. 

Pages on which the 35 Km layer can be seen are 93, 105, 103, 114, 117, 
126, 132, 141, 144, 150, 153 and T59. Only the scans on pages 153 and 159 
meet the criteria, and the fit for the data on page 159 (seen on page 157) 
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is suspicious because of the large deviation in the first half of the curve 
from the model and previous results. 

It seems that the apparent layers seen at 35 Km and 40 Km represent a 
single layer which is located in the inversion process at different altitudes 
due to inaccurate scaling. Because five of seven scans that placed the layer 
at 40 Km meet the above criteria and a sixth (No. 120) almost does but only 
two of twelve that placed it at 35 Km meet the criteria, we conclude that 
40 Km is most likely the correct altitude of the layer. 

We acquire more confidence in the 40 Km results when we note that in 
pairs the charts represent the same data scaled to the two models; that is, 
Charts 90 and 111, Charts 99 and 120, and Charts 135 and 156 are the results 
of scaling and inverting three data sets to the 1964 and 1968 models 
respectively. The first two sets are band 1 from pass 47 at two different 
locations (different scan numbers). The third data set is band 2 from 
pass 61 which, of course, is a third location. The numerical results between 
models remain significantly different. Within each model the results for the 
40- Km layer are consistent: the 1964 model places the attenuation coefficient 

ratio of aerosol to Rayleigh at .1; the 1964 places che value at .23. 

In reviewing the 20 Km layer one must account for the inherent error 
in the first few points caused by the inversion process. With this in mind, 
almost all of the SI 92 charts indicate a large 20 Km layer exists. Again, 
where we restrict our attention to those which satisfy the above criteria for 
a successful scale, charts 90, 99, 111, 135, and 156 emerge. And again we 
will include chart 120 which almost satisfies it. As before. Charts 90 and 
111 represent the results using models 1964 and 1968 respectively, as do 
Charts 97 and 120 and Charts 135 and 156, for the same measured data. Each 
of the pair of charts represent a different location (scan number or pass). 

The numerical results are model dependent but consistent within each model: 
the 1964 model finds the attenuation coefficient ratio (aerosol :Ray1eigh) 
to be about .1; the 1968 model finds it to be between .75 and 1.0. 

SI 91 Results 

The S191 was a much more sensitive Instrument than the SI 92 but 
apparently its absolute radiometric calibration was less accurate. Hence, 
the criteria for selecting the best SI 91 scaling and inversion results vyere 
changed from the SI 92 criteria as follows: 
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1) .75 < a < 1.25 

2) a good agreement in brightness values between the mode! and the 
measurements for the upper two thirds of the curve 

3) No defined restriction on the bias value, b, because a deep space 
value was not obtained. 

Where the 5192 instrument was sensitive in the lower wavelengths to about 45 
Kjti, the 5191 instrument was sensitive to about 70 Km. Hence, we could re- 
quire that a larger portion of the brightness curve match the model since 
the model is less influenced by aerosol changes in these regions. Also, the 
restrictions on the multiplicative scale factor, a, was relaxed as a result 
of the documentation stating that 5191 radiometric values could be in 
error by 15?i. The only restriction on the bias value, b, was that it be 
small in comparison to the data values. 

Aerosol layers are apparent in many bands at about 20, 40, 50, 60, and 
67 Km. Those at 20 and 40 Km are consistent with the 5192 results. Charts 
160-201 contain 5191 charts of the results. All show several layers con- 
sistent among themselves. Charts 162, 165, 174, and 183 are those which 
satisfy the above criteria. Table 6 summarizes the information on these 
charts showing for each chart the location of an aerosol layer and the value 
of the attenuation coefficient ratio (aerosol :Rayleigh) at that location. 

Within the 1964 model (No. 162, 165, and 174) is an indication of the 
wavelength dependence of the layer, and hence its aerosol size. However, 
one must be careful not to conclude too much from such few data points. 

It would be valuable to extract more points from the 5191 data in the future 
to-better derive conclusions. From the information in table 6 we note that 
layers at 40, 50, and 55 Km are more responsive to the longer wavelengths 
shown while layers at 59 and 66 Km are most responsive around .53 n falling 
off at higher and lower wavelengths. 

Charts 162 and 183 are the results of scaling and inverting the same 
data to the 1964 model and the 1968 model , respectively. Note that the 
same layers are evident but the derived attenuation coefficient values, 
and hence the ratios (aerosol ;Rayleigh) are different. For example, 
respectively for the 1964 model and the 1968 model ratio values obtained are 
.07 and .24 at 23 lOn, .04 and .22 at 39 Km, .09 and .26 at 48 Km, and .7 
and .86 at 59 Km. 
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Approximate Altitude 5 


Chart No. 


(Page No. ) 


15 

20 

30 

40 

50 

55 

60 

65 

ra 

rs 

o 

S 

162 

a) 


23 

29 

39 

48 

55 

59 

66 

.48 

64 


b) 


.07 

.05 

.04 

.09 

.16 

.7 

.04 



165 

a) 


19 


38 

49 


59 

67 

.53 

64 


b) 


-.2 


.2 

.4 


1. 

.7 



174 

a) 

14 

19 

29 

41 

47 

52 

59 

63 68 

,71 

64 


b): 

.1 

.08 

.04 

.24 

.58 

.43 

.43 

,2 .04 



183 

a) 


23 

30 

39 

48 

55 

59 

67 

.48 

68 


b) 


.24 

.2 

.22 

.26 

.37 

.86 

.3 




Table 6. SI 91 Summary Results 

(a) Aerosol Layer Altitude 

(b) Aerosol .‘Rayleigh attenuation coefficient ratio 

S19QA Results 

The results from the ST90 sensor must be presented as only qualitative 
since the attempts to scale and calibrate the data were unsatisfactory. This 
is undoubtedly because of the limited sensitivity of the duplicate film. 

The film was sensitive to light intensities below 30 Km where the scaling was 
severely model dependent. We therefore report the findings as relative and 
qualitative. The altitude and coefficient numbers are only for comparison 
between charts with no real v/orld meaning. The data was scaled against 
two models which are labeled 1964 and 1968 models. These are not the same 
models used for the SI 91 and SI 92 data analysis, however, because of a 
programming error. It was felt that there was little advantage In redoing 
the work since the results were not quantitative anyway and the models 
served well enough to assign an approximate altitude sufficient for in- 
version. 

For station 6, X = .55 y, pages 204, 207, and 210 show the results 
from the 1964 model. Layers are apparent at 16, 22, and 26 Km, On pages 
219, 222, and 225 the same data scaled to the 1 968 model shows layers at 
22, 26, and 32 Km. 

For station 5, X = ,68 y, pages 213 and 216 show the results from the 
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1964 model. Layers are apparent at 15, 20, and 22 Km. The 1968 model results 
are shown on pages 228, 231, and 234. Page 237 contains results from station 
2, A = .85 ,ini scaled to the 1968 model. 
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Table 7. Multiplicative and Additive (xlO"^) Scale Factors 
for si 92 Dataj Passes 61 and 47 


PASS 61 


1964 Model 


1968 Model 


Dark Space 
Bias Value 
-4< 


(100-124) (600-624) (120-144) n00--124) (600-624) (120-144) (x 10 


.9576 

1.267 


1.0573 .9202 .9507 

9.174 .6333 -1.242 


.9802 - 1.0726 

-1.200 - -2.377 


1.105 1.0294 

7.214 -1.716 


1 .0621 1.1227 1.1087 

4.049 -.9139 -1.690 


PASS 47 


1 


(100-124) 

(350-374) 

(100-124) 

(350-374) 


! 

1 

•1 

a 

.9990 

.9956 

.9545 

.9522 


b 

-.8209 

-.5310 

-.520 

-.3153 

.889 

1 - 

i 9 

a 


i.Ton 


1.0608 


c 

b 

- 

-3.284 

-■ 

-2,636 

1.583 

1 

■5 

a 

.9955 

1 .0700 

.9909 



: o 

b 

-3.475 

-3.104 

2.793 

- 

1.414 

> A 

a 



1.0255 



! 4 

b 

- 

■ - i 

-4.967 

- 

1.836 


a 

.9992 

.. 

_ 



i b 

J 

b 

-2.874 

- 

" 


,959 

- "7 

a 


1.0245 


T.07S2 


^ 7 

b 


-3.333 


-3.177 

1.406 


i 
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Table 8. Multiplicative Scale Factors for SI 91 Data, Pass 47 



1968 Model 

1964 Model 

.48 

1.17 

1,249 

.53 

.9283 

.9904 

.63 

.5006 (x) 

.5636 (x) 

.68 

.7477 

.8978 

.73 

.8281 

1.228 

.78 

.8607 

1.611 (x) 

.83 

.8264 

2.84 (x) 
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5.0 SUMMARY AND RECOMMENDATIONS 

The investigation results verify the analysis system developed to reduce 
limb aerosol data. The system is capable of analyzing large volumes of data 
quickly to yield altitude profiles of attenuation coefficients for each wave- 
length band. The sensors appear able to penetrate the atmosphere to about 
10 Km for the longer wavelengths and to about 20 Km for the shorter wave- 
lengths, The S192 mul tispectral scanner appears sensitive to an altitude 
of about 45 Km for the shorter wavelengths. The SI 91 sensitivity appears to 
extend to about 70 Km. The S190A duplicate film appears to be sensitive to 
about 29 Km. The data contained noise which had to be analyzed and removed 
which subtracted from the time available for aerosol analysis. In general, 
the sensors appeared adequate to detect stratospheric aerosols although 
future instruments would benefit from the experience gained through EREP. 

For example, more effort could be devoted to the pointing accuracy of the 
instruments thereby avoiding the necessity to compare with models in order 
to locate the field of view. Accurate radiometric calibration is also 
critical. Therefore efforts to attain the optimum calibration should be 
continued. 

The problems with the data and the effort it required to solve the 
problems prevented the analysis of the data for aerosol size distribution 
and number density in the manner described in the text. This would be a 
candidate for future work. 

In Section 4.0 much analysis results were presented. It was noted that 
several aerosol layers could be identified with a quantitative measurement 
of their attenuation coefficients. However, it was also noted that the 
results were dependent on the model invoked to scale and locate the data. 

The implication is clear that additional refinement of the models is neces- 
sary, also a candidate for future work. The point should be emphasized that 
the results were only as accurate as the model approximates the stratosphere. 
This could easily be the cause of inconsistencies. For example, the con- 
sistent increase of the values of attenuation coefficients in the lower 
altitudes above those predicted by the models is most likely due to the 
single scattering approximation when in fact multiple scattering is actually 
occurring in this region. Another improvement in the model would be to 
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include the effects of a scattering function (P{x)) (Section 2.1.1) that 
is dependent on altitude and particle size distribution. If the calibration 
of the data were inaccurate, a similar error would result in scaling and 
inverting the data. 

Much of the acquired data was not analyzed due to the lack of time 
caused by data complications, some of which were mentioned above. With the 
effort devoted to the development of a fairly involved system and with the 
understanding already gained about the data, continued analysis would be 
cost effective if stratospheric aerosols remain of interest. Additional 
effort to improve the models v^ould also be advisable since any improvements 
are directly extended to the results. 

The investigation has unquestionably verified that the described 
approach can contribute to the analysis of the particulate content of the 
stratosphere. The approach appears especially suitable to a future satellite 
which would monitor the changes and variations in the stratosphere between 
10 Km and 70 Km. The models used in much of this investigation would in 
that case become previous measurements. Differences would be interpreted 
as variations in time or space. 
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